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ABSTRACT 


Moment-curvature  relationships  have  been  developed  for  plates  with 
various  transverse  cross  sections  by  previous  investigators.  The  main 
purpose  of  this  thesis  was  to  extend  this  work  to  include  plates  with 
cambered  bitrapezoidal  cross  sections..  Experimental  results  are  presented 
showing  good  agreement  between  theory  and  experiment.  The  theory  indicat¬ 
ed  that  if  the  camber  was  increased  beyond  a  critical  value,  the  plates 
exhibited  a  bending  instability  for  certain  values  of  curvature.  This 
instability  characteristic  was  verified  by  experiment. 

The  theory  presented  was  for  a  general  cross  section.  By  changing 
certain  parameters  the  theory  was  shown  to  approach  the  linear  theory 
for  a  flat  plate  with  a  rectangular  cross  section  and  it  also  approached 
the  results  obtained  for  a  plate  with  a  bitrapezoidal  cross  section  with 


no  camber. 
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INTRODUCTION 

The  governing  equation  for  the  moment-curvature  relationship 
for  a  plate  with  a  rectangular  cross  section  can  be  written  as 

M  =  -  EI_  (1  -  ku^)  where  El  is  the  stiffness,  R  is  the  radius  of 
R  (1  -  y2) 

curvature  in  the  longitudinal  direction,  y  is  Poisson's  ratio  and  k 

is  a  parameter  which  depends  on  the  ratio  d^_  (d  is  the  width  and  t  is  the 

Rt 

thickness  of  the  plate).  For  a  beam  k  =  1,  and  for  a  plate  k  =  0  (i.e. 
a  plate  is  one  whose  width  is  many  times  larger  than  its  thickness). 

For  beams  and  plates  with  rectangular  cross  sections  there  is  a  gradual 
transition  in  k  from  a  value  of  one  to  zero  as  the  width  of  the  beam 
is  increased  while  the  thickness  is  held  constant. 

When  plates  with  tapered  cross  sections  are  considered,  the  moment- 
curvature  relationship  becomes  more  complicated .  The  aim  of  this  thesis 
was  to  investigate  theoretically  and  experimentally  the  stability  of 
plates  with  a  general  bitrapezoidal  cross  section  g.nd  to  show  that  with 
suitable  parameter  changes,  the  theory  approached  the  results  obtained 
by  other  investigators. 
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CHAPTER  I 

THEORETICAL  CONSIDERATIONS 

1.1  HISTORICAL.  REVIEW 

The  problem  of  finding  a  moment-curvature  relationship  for  a  flat 

plate  with  a  rectangular  cross  section,  with  its  longitudinal  edges 

o\* 

being  free,  was  first  investigated  by  Kelvin  and  Taitv  '  in  1879. 
Timoshenko*^)  modified  the  form  of  Kelvin  and  Tait's  equation  by 
including  the  k  term  which  accounted  for  the  transition  from  a  beam 
to  a  plate.  In  1949,  Fliigge  obtained  moment-curvature  relationships 
for  various  forms  of  tapered  cross  sections.  Some  of  the  sections  he 
considered  were  rectangular,  double-wedge  and  certain  bitrapezoidal 
cross  sections.  Fliigge  obtained  his  results  by  assuming  a  power  series 

(4) 

form  of  solution.  Ashwell  ,  in  1952,  considered  the  bending  instab¬ 
ility  of  curved  rectangular  plates.  Later,  Ashwell extended  his 

work  to  include  bending  instability  of  slightly  corrugated  rectangular 

( 

plates.  Fung  and  Wittricx  in  1954,  mentioned  the  problem  of  bending 
instability  for  plates  with  a  double-wedge  cross  section,  however,  no 
theory  was  presented. 

(3) 

One  of  the  cross  sections  considered  by  Fliigge  had  finite  edge 

(6) 

thickness  but  no  camber  (see  Fig.  1.1).  Fung  and  Wittrick  considered 
a  cross  section  with  infinitely  thin  edges  including  a  camber  as  shown 
in  Fig.  1.2.  In  the  bending  of  plates  with  cross  sections  as  shown  in 
Fig.  1.3A,  torsional  instability  develops.  Goodier^  discussed  torsional 
instability  in  thin-walled,  constant  thickness,  open  section  channels, 
however,  this  was  beyond  the  scope  of  this  thesis. 


*  Numbers  in  parenthesis  refer  to  references  at  the  end  of  thesis. 
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The  theory  developed  in  this  thesis  is  for  a  tapered  cross  section 
having  finite  edge  thickpess  and  camber  (Fig.  1.3)  which  does  not  exhibit 
torsional  instability.  By  a  change  in  sign  for  the  parameter  w0,  the 
theory  may  be  applied  to  a  plate  with  a  cross  section  as  shown  in  Fig.  1.3A. 
However,  torsional  instability  is  not  taken  into  account  in  the  analysis. 


FIG-  1-1  FIG.  1.2 

FLUGGE’S  CROSS  SECTION  FUNG  AND  WITTRICK'S  CROSS  SECTION 


FIG.  1.3 


FIG.  1.3A 


AUTHOR'S  CROSS  SECTION 


ALTERNATIVE  CROSS  SECTION 
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1.2  DEVELOPMENT  OF  THEORY 


The  governing  differential  equation  was  obtained  by  considering 
the  equilibrium  of  a  plate  element  as  shown  in  Fig.  1.4.  The  general 
cross  section  considered  is  shown  in  Fig.  1.5  and  the  manner  in  which 
the  plate  was  loaded  is  shown  in  Fig.  1.6. 

The  assumptions  used  in  the  development  of  the  theory  were  as  follows 

a)  The  behavior  of  the  material  under  loading  was  perfectly  elastic 
and  obeyed  Hooke's  Law. 

b)  The  material  was  homogeneous. 

c)  The  slope  of  transverse  deflection  curve  was  small  compared  to  unity. 

d)  The  weight  of  the  plate  was  zero. 

Applying  equilibrium  conditions  to  the  plate  element  (Fig.  1.4) 
and  using  geometrical  and  Hooke's  Law  relationships  the  following  diff¬ 
erential  equation  was  obtained 


dx‘ 


Ef 


12(1  -  y2) 


dzw 

dx2 


R-w 


+ 


Ewt 


(R-w) 


=  0 


1.2-1 


where  the  thickness  t  of  the  transverse  section  is 


t  =  2t0(l  -  x)  +  c 


1.2-2 


b 

for  any  positive  x  and  c  is  the  edge  thickness.  The  transverse  deflec¬ 
tion  w  can  be  assumed  to  be  small  compared  to  R,  hence 

R  -  w  i  R. 


Eqn.  1.2-1  can  now  be  written  as 


dx" 


Et; 


d2w 


+ 


12(1  -  yz)  l  dxz  R 


+ 


Ewt 

R2 


=  0 


1.2-3 


Lamb^)  obtained  a  solution  to  the  alpove  equation  for  a  plate  of  con¬ 


stant  thickness  t. 


■ 


' 
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FIG.  1.4  PLATE  ELEMENT 


.  FIG.  1.5  GENERAL  CROSS  SECTION 


dx 


FIG 


1.6  PLATE  ORIENTATION 
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Using  von  Karman's  large  deflection  equations,  Fung  and  Wittrick 

/  r  \ 

obtained  the  same  result  as  Lamb.  These  authors v  extended  their  work 

to  include  the  case  where:  (a)  the  cross  section  of  the  plate  had 

a  variable  thickness  t(x), 

(b)  initial  mid-line  of  the  cross  section 
was  displaced  from  the  centroidal  axis, 
and 

(c)  the  cross  section  had  no  edge  thickness, 
(i.e.  c  =  0) . 

Bellow  and  Faulkner ^ extended  Fung  and  Wittrick’ s  work  to  a  cross 
section  (see  Fig.  1.3A)  which  included  (a)  and  (b)  above  and  an  edge 
thickness,  i.e.  c  >  0. 

To  simplify  calculations  let 

w  =  wQ  +  C,  1.2-4 

where  w0  is  the  initial  distortion  of  the  mid-line  of  the  cross  section 
from  the  neutral  axis,  and  £  is  the  additional  deflection  due  to  the 
applied  moments.  By  considering  the  first  moment  of  area  of  the  cross 
section  (Fig.  1.5)  about  the  neutral  axis,  it  can  be  shown  that 


w 


o 


t0(l  +  k) 


x_  -  3  +  2/3 
b  23  +  1 


where  ,n  c  . 


Since  w0  is  a  linear  function  of  x,  Eqn.  1.2-3  becomes 


1.2-5 


d±_  / 1 3  dh  j  +  V_  d2  (t3)  +  12t(l  -  li2)  (4  +  w0)  =  0 
dx  \  dx2/  R  dx2  r2 


1.2-6 


This  can  be  transformed  to  a  simpler  form  by  setting 


p  =  t 


o 


5  =  1- 


X 


and 


7 


yielding 


df_  fp  3  dfc  \+  d2  (p3) 
d£2  V  d£2  /  R  d£2 


+  4X40 (4  Hh  WD)  =  0 


1.2-7 


where  X4  =  jb-.!1  ~  d2)  and  t,_  is  defined  as  one-half  the  thiek- 

R2  tQ2 

ness  of  the  tapered  portion  of  the  plate  (Fig.  1.5).  Further  simplifica¬ 
tions  were  made  by  substituting 

Y  =  5  +  3, 

giving  the  final  form  of  Eqn.  1.2-7  as 


+  6yb2y  +  X4y^  +  wo^  =  °* 
R 


1.2-8 


In  the  solution  of  Eqn.  1.2^8  let  40  be  the  particular  solution 
and  4*  be  the  complimentary  solution.  Thus,  40  becomes 

s  =  -  „  -  6ubi,  1.2-9 

0  a4r 

and  the  complimentary  equation  becomes 

ll_  (y3  d2^V  X4y4*  =  0.  1.2-10 

dy2  \  dy2  ) 

Eqn.  1.2-10  is  of  the  same  form  as  obtained  for  the  solution  of  a 

(11) 

cylindrical  tank  with  non-uniform  wall  thickness.  Timoshenko  and 

(12) 

Fliigge  obtained  the  following  solution  to  Eqn.  1.2-10: 

£*  =  .L  [A  ber^n  4-  B  bei^q  +  C  ker^n  +  D  kei^n]  1.2-11 

where  ri  is  defined  as 

n  =  2Xy^. 

The  complete  solution  may  be  written  as 

4  =  40  +  4* 
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or  after  substituting  for  and  C* ,  C  becomes 


C  -  —  [A  ber"n  +  B  bei^n  +  C  ker^p  +  D  kei'n] 

n 


w  -  ^ b 2  1.2-12 

X4R 


where  the  prime  denotes  differentiation  with  respect  to 
The  constants  of  integration  A,  B,  C,  and  D  can  be 
following  four  conditions: 


n. 

solved  from  the 


(a)  ~  — 7i  ^  .,2- 


12(1  -  uO 


d2w  +  y_ 
dx2  R 


=  o, 


x=b 


(b)  Vx  = 


since 


Qx  +  9_(Mxy) 
3y 

Wx=b=  0> 


=  0 

x=b 

(b)  reduces  to 


rd3w" 


=  0  . 


_dx3Jx=b 


(c)  from  the  symmetry  of  the  problem 


(d) 


dx  =  0. 


and 


The  above  equations  were  transformed  into  non-dimensional  form  and 
the  constants  A,  B,  C,  and  D  were  solved  using  an  IBM  7040/1401  digital 
computer . 


1.3  MOMENT-CURVATURE  RELATIONSHIP: 

By  considering  moment  equilibrium  of  an  elemental  strip  of  width  dx, 
the  external  moment  dM  applied  to  the  transverse  edge  must  balance  the 
two  internal  moments  (My)dx  and  (Nyw)dx  (Fig.  1.7).  Writing  the  equili¬ 
brium  equation  for  Fig.  1.7 

dM  =  (My) dx  +  (Nyw)dx 


or  M  =  2  \  (My) dx  +  2  \  (Nyw)dx 


1.3-1 


noting  that  the  cross  section  is  symmetrical  about  the  centre  line.  The 
internal  moment  Nyw,  is  the  moment  per  unit  width  which  is  induced  when 
the  centre  line  of  the  cross  section  is  displaced  from  the  centroidal 
axis  by  an  amount  w. 


FIG.  1.7  PLATE  STRIP 


For  ease  in  evaluating  Eqn.  1.3-1*  it  is  advantageous  to  integrate 
with  respect  to  r\ . 

Since  the  slope  of  the  transverse  deflection  curve  is  small  compared 

to  unity,  the  curvature  in  the  transverse  direction  can  be  written  as 

2. 

1  =  ^  w  >  Eqn.  1.2-5  shows  that  wG  is  a  linear  function  of  x  and 

Rx  dx2 

w  =  w0  +  C  (Eqn.  1.2-4)  thus,  1/RX  can  be  written  as 


1_  =  d2C  . 

Rx  dx2 


1.3-2 
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The  equation  for  becomes 


My  =  ~Et 


12(1  -  y2) 


-  +  y 


d2C 


R  dx" 


where  1_  ^  1  . 
R  R 


1.3-3 


Since  ^  is  a  function  of  p,  then  by  applying  the  rules  of  differentiation 
the  following  is  obtained: 

d2 4  =  d2£  r dpi2  +  d2n  dc  .  1.3-4 

dx2  dp2  Ldx J  dx2  dp 


Let  the  equation  for  C  be  written  as 


C 


wr 


6yb2 

X4R 


1.3-5 


where  F'  is  defined  as 


F^  =  A  ber^p  +  B  bei^p  +  C  ker"p  +  D  kei^p  . 


1.3-6 


By  making  appropriate  substitutions,  Eqn.  1.3-4  becomes 


dh  = 

dx2 


—f'  -  —f"  +  If' 


t0(i  +  k) 

2A2 


^bp 


±-F'  -  — F 


+ 


ptp(l  +  k) 

2A2 


4A" 


1.3-7 


.  b  2  p ' V 


The  double  and  triple  prime  of  F  are  the  second  and  third  derivatives 
respectively  of  F  with  respect  to  p. 

The  membrane  force  Ny  can  be  calculated  from  the  membrane  stress> 


i>-  V7  i 
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yielding  t 

7 


N  =  \  a,,  dz' 

y  \-t 7 


1.3-8 


where  z'  is  measured  from  the  center  line  pf  the  cross  section  and  t  is  the 
thickness  of  the  plate  at  the  point  under  consideration.  Evaluation  of 
Eqn.  1.3-8  yields 


m  _  -Etw 
7  R~ 


but,  since  Ry  ±  R  then 


N  =  "Etw 
y  R  ' 


1.3-9 


From  Eqns .  1.3-1,  1.3-3,  1.3-7,  and  1.3-9,  the  following  was  obtained: 

'£ 


M 


Ebt0 


J 


96A8  (1  -  y: 

l)\ 

R  .  b 

/e 

'a 

+  bEt0  (  3 

v 

2A4R  \ 

n 

A4R 

1 

dn 


dp 


1.3-10 


The  limits  of  integration  e  and  a  are  obtained  from  the  relationship 

n  =  2 A ( 1  -  X  +  3)^  . 

b 

The  original  limits  of  integration  were  x  =  b  and  x  =  0.  Due  to  the 
transformation  this  yields  n  =  2X8^  and  n  =  2 A ( 1  +  3)  2  respectively. 
For  simplicity  e  and  a  are  defined  as 
e  =  2A (6)^ 

and  a  =  2X(1  +  8)^  . 


Integration^13^  of  Eqn.  1.3-10  leads  to  the  following: 


■ 
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M 


Ebt03 
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where  the  subscripts  on  the  preceedlng  Bessel  functions  indicate  the 
order  of  the  function.  No  subscript  indicates  a  Bessel  function  of 
order  zero.  The  solution  to  Eqn.  1,3-11  was  obtained  using  an  IBM 
7040/1401  digital  computer.  Details  of  this  solution  along  with  the 
solution  for  the  constants  of  integration  appear  in  Appendix  A. 

1.4  DISCUSSION  OF  THEORETICAL  RESULTS 

Theoretical  curves  for  the  moment-curvature  relationship  were 
obtained  for  different  parameter  variations.  The  effect  on  the  moment- 
curvature  relationship  when  k  was  increased  while  other  parameters 
were  held  constant  is  shown  in  Fig.  1,8  when  the  plate  was  turned  upside 
down  (concave  side  down) .  As  k  was  increased,  the  slope  of  the  curves 
increased  sharply  due  to  the  increased  plate  stiffness,  however,  a 
maximum  moment  was  reached  after  which  the  moment  decreased  for  an  in¬ 
crease  in  curvature.  The  point  of  instability  was  reached  when  the 
effect  of  the  membrane  and  transverse  bending  forces  were  large  enough 
to  flatten  out  the  cross  section  of  the  plate.  As  the  camber  decreased 
due  to  the  increased  loads,  the  plate  began  to  act  as  a  flat  plate. 
Theoretically,  for  large  curvatures,  the  moment-curvature  relationships 
for  different  k  values  should  approach  the  case  where  k  =  0,  This  is 
shown  to  be  the  result  in  Fig,  1.8. 

In  Fig.  1.9,  moment-curvature  relationships  were  plotted  for 
the  case  when  the  plate  was  bent  with  the  concave  side  up.  In  the 
cross  section  used  in  Fig.  1.8,  the  membrane  forces  and  the  trans¬ 
verse  bending  action  are  both  acting  to  distort  the  section  anti- 
clastically .  When  the  plate  had  the  concave  side  up  (Fig.  1.9),  the 
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FIG,  1.8  THEORETICAL  MOMENT^- CURVATURE  RELATIONSHIPS  FOR  VARIATION  IN 

(CONCAVE  SIDE  DOWN) 
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FIG.  1.9  THEORETICAL  MOMENT- CURVATURE  RELATIONSHIPS  FOR  VARIATION  IN 

(CONCAVE  SIDE  UP) 
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FIG.  1.10  THEORETICAL  MOMENT- CURVATURE  RELATIONSHIPS  FOR  VARIATION  IN 

(CONCAVE  SIDE  DOWN) 
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FIG.  1.11  THEORETICAL  MOMENT -CURVATURE  RELATIONSHIPS  FOR  VARIATION  IN 

(CONCAVE  SIDE  UP) 
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transverse  bending  action  distorted  the  plate  anticlastically  but  the 
components  of  the  membrane  forces  counteracted  this.  A  critical  point 
is  reached  when  the  effect  of  the  membrane  forces  becomes  larger  than 
the  effect  of  the  transverse  bending  action  and  the  cross  section  j 
flattens  out.  It  is  due  to  these  opposite  effects  that  the  maximum 
moment  reached  for  this  section  (Fig.  1.9)  is  greater  than  for  Fig.  1.8 
before  the  cross  section  distorts  into  a  flattened  configuration. 

Figs.  1.10  and  1.11  illustrate  the  effect  on  the  moment-curvature 
relationship  when  the  edge  thickness,  c,  is  varied  while  the  other 
parameters  are  held  constant.  As  expected,  the  value  of  the  moment 
increases  with  an  increase  in  c  for  a  particular  value  of  curvature, 
although  the  general  shape  of  the  curves  remains  the  same. 

1.5  COMPARISON  OF  THEORY  TO  PREVIOUS  RESULTS 

Flligge^  obtained  the  moment-curvature  relationship  for  a  plate 

(2) 

with  a  cross  section  of  the  type  shown  in  Fig.  1.1  and  Timoshenko v  y 
(14) 

and  Searle  obtained  the  moment-curvature  relationship  for  a  flat 
plate  with  a  rectangular  cross  section.  By  changing  certain  parameters 
in  the  theory  developed  in  Secs.  1.2  and  1.3,  both  of  the  cross  sections 
considered  by  the  authors ^ ^ were  obtained.  Fung  and  Wittrick^ 
stated  the  moment-curvature  relationship  for  a  plate  with  a  cross  section 
as  shown  in  Fig.  1.2,  however,  they  did  not  present  any  results  in  their 
paper. 

In  comparing  the  theoretical  solution  with  Flligge's  results,  the 
cross  section  had  a  finite  edge  thickness,  c  >  0,  and  a  finite  taper 
t0  >  0.  The  cross  section  and  the  results  are  shown  In  Fig.  1.12. 


19 


The  theoretical  solution  was  obtained  by  setting  k  =  0.  Fliigge '  s  results 
and  linear  theory  shown  are  only  approximate  as  they  were  obtained  from  a 
graph. 

In  comparing  the  developed  theory  with  the  theory  for  a  plate  with 
a  rectangular  cross  section,  t0  and  k  have  to  be  set  equal  to  zero.  As 
tQ  approaches  zero,  X  approaches  infinity  and  no  solution  can  be  obtained. 
For  comparison,  tQ  was  decreased  to  about  10  percent  of  the  edge  thickness. 
The  cross  section  considered  and  the  comparison  is  shown  in  Fig.  1.13. 

It  is  evident  that  if  t0  approached  zero  the  two  theories  would  yield 
the  same  results. 

Figs.  1.12  and  1.13  indicate  that  agreement  between  different  solu¬ 
tions  is  good.  In  the  following  sections,  experimental  verification  is 
presented  for  plates  with  different  variations  of  the  cross  section  as 
shown  in  Fig.  1.5. 


Moment,  In. -Lbs 


20 


FIG.  1.12 


COMPARISON  OF  THEORY  WITH  PREVIOUS  WORK  FOR  BITRAPEZOIDAL  CROSS  SECTION 
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FIG.  1,13  COMPARISON  WITH  LINEAR  THEORY  FOR  A  RECTANGULAR  CROSS  SECTION  PLATE 
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CHAPTER  II 

EXPERIMENTAL  APPARATUS 


2.1  PLATES 

The  plates  tested  were  machined  from  a  flat  rolled  sheet  of  Alcoa 
7075-T6*  Aluminum.  The  modulus  of  elasticity  of  this  material  was  10.4  x 
10  psi  with  a  Poisson's  ratio  of  0.333.  With  a  lower  modulus  of  elas¬ 
ticity  than  for  steel,  higher  curvatures  were  obtained  with  smaller  loads 
than  would  be  possible  using  steel.  The  density  of  aluminum  being  lower 
than  for  steel,  allowed  an  aluminum  plate  to  act  more  as  a  weightless 
plate  than  a  steel  plate  of  equal  dimensions. 

Plate  I-A  was  machined  from  a  sheet  of  aluminum  with  a  nominal  size 

of  0.125  x  12.0  x  67.0  inches.  As  shown  in  Fig.  2.1,  the  cross  section 

warped  during  machining.  For  calculation  purposes,  the  cross  section 

was  assumed  to  have  a  linear  thickness  variation.  The  final  cross  section 

could  be  bent  to  an  approximate  longitudinal  curvature  of  0.100  inches  ^ 
without  yielding  the  material. 

After  the  necessary  experiments  were  performed  on  Plate  I-A,  Plate 
I-B  was  obtained  by  performing  a  "controlled  bend"  along  the  length  of 
Plate  I-A.  The  bend  was  obtained  by  placing  Plate  I-A  in  a  press  brake 
and  then  applying  a  line  load  along  the  centre  line  throughout  the  length 
of  the  plate.  Measurements  were  made  after  every  load  increment  to 
ensure  that  the  plate  bent  uniformly  throughout  its  length.  The  final 
cross  section  is  shown  in  Fig.  2.2. 

^Properties  of  alloy  used  were  obtained  from  Alcoa  Structural  Hand  Book 
Aluminum  Company  of  America  1956,  p.  35. 
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Axis  of  Symmetry 
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FIG.  2.1  CROSS  SECTION  OF  PLATE  I-A 
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FIG.  2.2  CROSS  SECTION  OF  PLATE  I-B 
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FIG.  2.3  CROSS  SECTION  OF  PLATE  I-C 
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Plate  I-C  was  obtained  by  turning  Plate  I-B  upside  down.  The  cross 
section  is  shown  in  Fig.  2.3. 

Plate  II-A  was  also  machined  from  a  sheet  of  7075-T6  aluminum.  The 
section  warped  slightly  during  machining;  however,  as  seen  in  Fig.  2.4, 
the  warping  was  quite  linear  and  an  accurate  description  of  the  actual 
cross  section  could  be  made..  The  final  cross  section  could  be  bent  to 
a  longitudinal  curvature  of  approximately  0.054  inches  ^  before  yeild- 
ing  would  occur. 

Plate  II-B  was  obtained  in  the  same  manner  as  Plate  I-B  (i.e.  by 
performing  a"controlled  bend"  along. the  length  of  the  plate). 

By  turning  Plate  II-B  upside  down,  Plate  II-C  was  obtained.  The  dif¬ 
ference  in  some  of  the  dimensions  between  Plates  II-B  and  II-C  are 
attributed  to  the  cross  section  of  Plate  II-B  flattening  out  during 
the  process  of  performing  experiments  on  the  Plate.  Hence,  the  cross 
section  for  Plate  II-C  had  less  camber  than  Plate  II-B.  The  cross 
sections  are  shown  in  Fig.  2.5  and  Fig.  2.6  for  Plates  II-B  and  II-C 
respectively. 

2.2  STRAIN  GAUGES 

Experimental  longitudinal  curvatures  were  calculated  from  the  bend¬ 
ing  strains  induced  in  the  top  and  bottom  surface  fibres  of  the  plate. 

These  strains  were  measured  with  Budd  Instruments  type  C12-121  strain 
gauges.  The  C12-121  is  an  electrical  resistance  metal  film  type  gauge. 

These  gauges  are  temperature  compensated  and  have  a  gauge  length  of  0.122 
inches  and  width  of  0.125  inches.  The  gauges  were  mounted  on  the  centre  line' 
at  various  positions  along  the  longitudinal  axes  of  the  plates.  (See 
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Axis  of  Symmetry 
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.  .  r  .  _  FIG.  2.5  CROSS  SECTION  OF  PLATE  II-B 
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FIG.  2.6  CROSS  SECTION  OF  PLATE  II-C 
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Fig.  2.7,  2.8,  2.9,  2.10).  The  mounting  was  such  that  the *top  gauges 
were  directly  above  the  ones  mounted  on  the  bottom  of  the  plate. 

Two  main  advantages  of  these  gauges  were  the  low  transverse 
sensitivity  and  the  thin  backing  material.  Since  the  longitudinal 
strains  were  of  interest  only,  then  by  using  low  transverse  sensitivity 
gauges,  the  transverse  strains  had  very  little  effect  on  the  longitu¬ 
dinal  strain  readings.  Thinner  backing  material  allowed  the  gauge  to 
be  placed  closer  to  the  surface  for  which  a  reading  is  desired,  and 
hence,  a  more  accurate  surface  strain  could  be  obtained. 

2.3  LOADING  APPARATUS 

A  uniformly  distributed  moment  was  applied  to  the  transverse  ends 
of  the  plate  by  means  of  bolting  the  two  ends  of  the  plate  to  two 
separate  one-inch  diameter  shafts.  A  fifteen- inch  diameter  pulley  was 
fastened  to  each  end  of  the  shaft.  By  means  of  wire  cables,  hangers 
were  hung  from  the  rims  of  the  pulleys  and  the  desired  loads  were  placed 
on  the  hangers. 

To  reduce  the  frictional  effects  at  the  supports,  the  shaft- 
pulley  system  was  placed  on  bearings.  The  bearings  rested  on  two  channels 
whose  surfaces  were  machined  smooth.  To  further  reduce  the  effects  of 
friction,  after  each  load  increment  was  applied,  the  shaft-bearing 
system  was  given  a  small  horizontal  displacement  on  the  supporting  chan¬ 
nels  and  then  it  was  allowed  to  come  to  equilibrium  before  any  strain 
readings  were  recorded.  For  a  photographic  view  of  the  apparatus,  see 


Fig.  2.11. 
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1  Measuring  Station 


Cross  Section  of  Plate  I-A 


FIG.  2.7  MEASURING  STATION  FOR  PLATE  I-A 


Cross  Section  of  Plate  I-C 


FIG.  2.8  MEASURING  STATIONS  FOR  PLATES  I-B  AND  I-C 


28 


Mid-Line 


T~ 

f 

K-10" - 

■  Jr5"  i 

i  n  i*a 

—  *  a  1 

-5'J1 
- 10"- — > 

-* - 38 

" - > 

■  Measuring  Station 


Cross  Section  of  Plate  II-A 


FIG.  2.9  MEASURING  STATIONS  FOR  PLATE  II-A 


Cross  Section  of  Plate  II-B 


FIG.  2.10  MEASURING  STATIONS  FOR  PLATES  II-B  AND  II-C 
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FIG.  2.11  PHOTOGRAPHIC  VIEW  OF  LOADING  APPARATUS 
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2  .  4  MEASUREMENT  OF  CURVATURE 


The  longitudinal  curvatures  were  calculated  from  the  strains  in¬ 
duced  in  the  top  and  bottom  fibres  of  the  plate.  A  BLH  switching  and 
balancing  unit  was  used  to  switch  from  one  measuring  station  to  the 
next  and  a  BLH  strain  indicator  was  used  to  read  the  strains. 

For  a  plate  in  pure  bending,  it  can  be  shown  that 

I  =  ~2 (ei} B  2.4-1 

R  t± 

where  (e^)g  =  bending  strain  at  x-j_ 
and  t^  =  thickness  of  the  plate  at  x^. 

However,  when  a  plate  is  bent  into  a  longitudinal  curvature  such  that 
axial  strains  are  developed  in  addition  to  bending  strains,  it  becomes 
necessary  to  eliminate  the  axial  strains  before  the  curvatures  can  be 
calculated.  The  axial  strains  were  eliminated  by  having  the  gauges 
wired  in  series  such  that  the  strains  at  the  top  of  the  plate  and  the 
strains  directly  at  the  bottom  of  the  plate  were  additive  giving  a 
double  output  of  the  actual  strain  on  the  top  or  bottom  surface. 
Defining  the  double  output  of  the  bending  strain  as  eB  and  substitut¬ 
ing  this  value  for  2(e^)g  in  Eqn.  2.4-1  the  curvature  becomes 

1  -eB  2.4-2 


Using  Egn.  2.4-2,  a  computer  program  was  written  to  evaluate 
curvatures  from  the  strain  gauge  readings  (see  Appendix  B) .  Knowing 
both  the  moment  and  the  corresponding  curvature,  the  experimental 


curves  were  plotted. 
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2 . 5  TESTING  PROCEDURE 


The  testing  procedure  used  for  all  six  plates  was  essentially  the 
same.  Initially,  a  flat  support  was  inserted  under  the  plate.  With 
the  flat  support  in  place,  the  strain  gauges  were  zeroed.  Then  the 
flat  support  was  removed,  allowing  the  plate  to  sag  due  to  its  own  weight. 
Loads  were  then  applied  in  desired  increments.  The  applied  load  and 
the  strain  readings  were  recorded  for  each  measuring  station.  After 
the  maximum  load  had  been  applied,  all  the  load  was  removed  and  the 
flat  support  was  reinserted  under  the  plate  so  that  the  zero  readings 
could  be  checked. 
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CHAPTER  III 


EXPERIMENTAL  AND  THEORETICAL  RESULTS 


3.1  REMARKS  ON  EXPERIMENTAL  RESULTS 

Three  different  sets  of  strain  gauge  readings  were  obtained  for 
each  measuring  station.  The  variation  in  strain  for  any  given  curvature 
was  small,  it  was  decided  that  an  average  obtained  from  three  sets  of 
readings  was  sufficient  to  establish  the  experimental  curve.  Each  experi¬ 
mental  plate  was  symmetrical  about  its  longitudinal  axis  and  the  measuring 
stations  were  placed  along  this  axis.  Since  the  minimum  number  of  stations 
on  any  plate  was  three  and  the  maximum  was  six,  this  gave  a  minimum 
number  of  nine  and  a  maximum  of  eighteen  different  readings  being  averaged 
for  each  point  plotted. 

3.2  COMPARISON  BETWEEN  EXPERIMENTAL  AND  THEORETICAL  RESULTS 

In  the  development  of  the  theory,  it  was  assumed  that  the  plate  was 
weightless;  hence,  when  the  experimental  and  theoretical  moments  were 
compared  for  a  particular  longitudinal  curvature,  it  was  necessary  to 
subtract  from  the  applied  moment  the  moment  due  to  the  weight  of  the 
plate.  It  can  be  seen  from  Fig.  3.1  that  the  weight  per  unit  horizontal 
distance  varies  with  the  distance  from  the  mid-line  at  a  particular  long¬ 
itudinal  curvature.  This  variation  must  be  accounted  for  in  the 
correction  moment.  By  summing  moments  about  a  section  p  units 
from  the  mid-line,  the  correction  moment  due  to  the  weight  of  the  plate 
can  be  shown  to  be 


3.2-1 
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where  M^  =  correction  moment, 

q  =  weight  per  unit  length  of  plate, 
p  =  horizontal  distance  from  mid-line  to  section 
where  correction  moment  is  desired, 

R  =  longitudinal  radius  of  curvature  (average  R) , 
L  =  one-half  the  total  length  of  plate,  and 
a  =  one-half  the  distance  between  supports. 


Mid-iLine 


FIG.  3.1  DIAGRAM  FOR  OBTAINING  CORRECTION  MOMENT 
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The  experimental  curves  were  plotted  and  compared  with  theory  in 
Fig^.  3.2  to  3.7. 

In  measuring  k  for  the  different  plates,  it  was  found  that  this 
value  varied  by  about  10  percent  between  different  points  on  the  same 
plate.  In  calculating  the  theoretical  curves,  an  average  value  for  k 
was  used.  It  was  found  that  a  small  change  in  k,  for  k  <  1,  pro¬ 
duced  a  negligible  change  in  the  theoretical  curves  but  a  small  change 
in  k,  for  k  >  5,  produced  large  changes  in  the  curves.  Not  being 
able  to  measure  k  accurately  for  the  plates  with  k  >  5  was  a  pos¬ 
sible  source  of  error. 

3.3  DISCUSSION  OF  RESULTS 

As  seen  from  Figs.  3.2  to  3.7,  better  agreement  was  reached  be¬ 
tween  theory  and  experiment  for  the  tests  involving  the  Plate  II  series 
than  was  obtained  from  the  tests  for  the  Plate  I  series.  Better 
agreement  was  reached  because  during  the  machining  of  the  cross  section, 
the  Plate  II  series  warped  linearly  whereas  for  the  Plate  I  series, the 
warping  was  non  linear  (see  Figs.  2.1  to  2.6).  Hence,  a  better  ap¬ 
proximation  to  the  actual  cross  section  could  be  made  for  the  Plate  II 
series . 

Fig.  3.2  shows  the  experimental  and  theoretical  curves  for  Plate 

1  -  A.  A  24  percent  discrepancy  between  the  two  curves  is  shown  to 

occur  at  a  curvature  of  0.015  inches  \  but  the  error  reduced  to  only 

-1 

2  percent  at  a  curvature  of  0.046  inches  .  One  possible  reason  for 
the  better  agreement  at  higher  curvatures  was  due  to  the  distortion 
of  the  cross  section  at  higher  loads.  It  was  assumed  that  as  the 
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FIG.  3.2  COMPARISON  OF  THEORETICAL  AND  EXPERIMENTAL  RESULTS  FOR  PLATE  I-A 


Theoretical  Curve 
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FIG.  3.3  COMPARISON  OF  THEORETICAL  AND  EXPERIMENTAL  RESULTS  FOR  PLATE 
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FIG.  3.4  COMPARISON  OF  THEORETICAL  AND  EXPERIMENTAL  RESULTS  FOR  PLATE 


Experimental  Points 
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PIG  3.5  COMPARISON  OF  THEORETICAL  AND  EXPERIMENTAL  RESULTS  FOR  PLATE  II- 


Experimental  Points 
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FIG.  3.6  COMPARISON  OF  THEORETICAL  AND  EXPERIMENTAL  RESULTS  FOR  PLATE  II- 
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FIG  3.7  COMPARISON  OF  THEORETICAL  AND  EXPERIMENTAL  RESULTS  FOR  PLATE  II- 
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cross  section  distorted,  it  approached  the  cross  section  used  to  ob¬ 
tain  the  theoretical  curve,  giving  better  agreement  between  theory  and 
experiment  at  higher  curvatures. 

The  curves  for  Plate  I  -  B  are  shown  in  Fig.  3.3  A  38  percent  dis- 

-1 

crepancy  occured  in  the  region  where  the  curvature  was  0.004  inches 
This  discrepancy  was  attributed  to  the  chacteristic  torsional  instabil¬ 
ity  which  occurred  in  the  bending  of  plates  with  cross  sections  as 
shown  in  Fig.  3.3.  When  the  cross  section  of  the  plate  flattened  out 
sufficiently  so  that  torsional  bending  no  longer  dominated,  reasonably 
good  agreement  between  theory  and  experiment  was  reached.  At  a 
curvature  of  0.029  inches  \  the  error  was  14  percent  and  at  a  curvature 
of  0.049  inches  \  the  error  was  only  7  percent. 

Bending  a  plate  with  a  cross  section  as  shown  in  Fig.  3.4  did  not 

exhibit  any  torsional  instability.  Consequently,  good  agreement  between 

theory  and  experiment  was  reached.  According  to  theory,  the  plate  be- 

comes  unstable  when  the  curvature  is  approximately  0.0055  inches 
The  difference  between  theory  and  experiment  in  this  region  was  only 

6  percent.  When  the  next  increment  of  moment  was  applied  to  the 
plate,  the  cross  section  of  the  pla,te  flattened  out  instantaneously, 
thus  exhibiting  a  characteristic  condition  of  instability.  The 
flattening  of  the  cross  section  caused  the  curvature  to  increase 
from  0.0055  inches  ^  to  0.0285  inches  ^  before  this  additional  in¬ 
crement  of  moment  could  be  supported.  After  passing  through  this 
region  of  instability,  the  plate  behaved  more  linearly.  At  a  cur¬ 
vature  of  0.050  inches  ^  the  agreement  between  theory  and  experiment 
was  within  5  percent.  It  is  interesting  to  note  that  Ashwell^^  verified 
his  theory  by  carrying  out  experiments  on  steel  measuring  tapes,  obtain¬ 
ing  an  error  from  4  to  19  percent  below  that  predicted  by  the  theory. 
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Better  agreement  was  reached  for  the  Plate  II  series  than  for  the 
Plate  I  series.  Plate  II  -  A  {Fig.  3.5)  showed  an  increase  in  the  de¬ 
viation  between  theoretical  and  experimental  curves  with  increased  curv¬ 
ature;  however,  the  difference  was  only  8  percent  at  a  curvature  of 
-.1 

0,036  inches  .  Plate  II  -  B,  Fig.  3.6,  exhibited  the  same  form  of 
torsional  instability  as  Plate  I  -  B. 

Fig.  3.7  for  Plate  II  -  C,  showed  excellent  agreement  between  theory 
and  experiment.  This  plate  did  not  exhibit  the  bending  instability 
characteristic  shown  by  Plate  I  -  C.  Before  this  plate  could  exhibit 
bending  instability,  an  increase  in  k  would  have  to  be  made. 

3.4  CONCLUDING  REMARKS  ON  CHAPTER  III 

It  was  observed  that  for  all  six  plates,  the  experimental  moments 
were  less  than  or  equal  to  the  theoretical  moments  at  a  particular  curv¬ 
ature.  One  possible  explanation  for  this  was  that  the  theoretical  curves 
were  based  on  the  dimensions  of  the  plate  before  the  plate  was  bent  into 
a  longitudinal  curvature.  During  the  process  of  conducting  three  differ¬ 
ent  tests  on  each  plate,  the  plates  gradually  lost  their  camber;  hence, 
reducing  their  stiffness  and  decreasing  the  moment  required  to  bend  them. 

In  addition  to  the  errors  caused  by  the  non-linearity  warping  of 
the  cross  section  and  the  inability  to  measure  k  accurately,  other  sources 
of  error  are: 

(1)  error  in  reading  the  strain  gauges,  and 

(2)  errors  due  to  the  frictional  effects  at  the  supports 
of  the  loading  apparatus. 

The  error  in  reading  the  strain  values  was  less  than  *  5  percent, 
and  the  frictional  error  was  approximately  i;  2  percent.  The  total 
error  due  to  the  above  causes  was  less  than  *7  percent. 
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CHAPTER  IV 

CONCLUSIONS  AND  FURTHER  PROBLEMS 

4.1  CONCLUSIONS 

Theoretical  and  experimental  investigations  were  described  in  this 
thesis  for  moment-curvature  relationships  of  cambered  bitrapezoidal 
plates.  The  theory  presented  was  supported  by  experimental  results. 

The  general  theory  developed  was  shown  to  approach  the  previous  theories 
for  flat  rectangular  plates  and  for  bitrapezoidal  plates  with  no  camber. 
In  fact,  the  previous  theories  are  special  cases  of  the  theory  presented 
in  this  thesis. 

In  view  of  the  good  agreement  obtained  between  theory  and  experiment 
and  with  previous  theories  for  different  cross  sections,  the  theory  can 
be  expected  to  accurately  predict  the  moment-curvature  relationships  for 
plates  with  cambered  bitrapezoidal  cross  sections  having  an  edge  thick¬ 
ness  . 

4.2  PROBLEMS  FOR  FURTHER  STUDY 

In  conducting  experiments  on  the  plates,  it  was  observed  that  Plates 
I-B  and  II-B  exhibited  torsional  instability.  It  would  be  interesting 
to  investigate  experimentally  and  theoretically  the  torsional  instability 
of  cambered  bitrapezoidal  plates. 

Thus  far,  plates  have  been  considered  with  linear  thickness  varia¬ 
tions"  The  present  theory  could  be  extended  to  include  thickness  varia¬ 
tions  of  second  or  higher  order  (i.e.  semi-elliptic  or  semi-circular 
cross  sections  with  camber) . 
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APPENDIX  A 

COMPUTER  PROGRAM  TO  EVALUATE  CONSTANTS  OF  INTEGRATION  AND  THEORETICAL 

MOMENT-CURVATURE  RELATIONSHIPS 


The  following  program  consists  of  three  parts:  main  program 
called  "PLATE"  and  two  subroutines  called  "SOLVE"  and  "BSL" .  First, 
the  main  program  was  used  to  solve  for  the  order  zero  and  one  Bessel 
functions  ber,  bei,  ker,  and  kei.  By  using  various  combinations  of 
the  order  zero  and  one  Bessel  functions,  the  order  two  and  the  deriva¬ 
tives  of  the  Bessel  functions  were  calculated.  A  recursion  formula 
was  used  to  evaluate  the  above  order  zero  and  one  functions,  however, 
for  arguments  greater  than  ten,  it  was.  discovered  that  this  recursion 
formula  gave  erroneous  results  for  the  ker  and  kei  functions  of  order 
zero  and  one.  Then  subroutine  "BSL"  was  added  to  the  program  so  that 
whenever  the  argument,  n> 10 ,  this  subroutine  was  called  upon  by  the 
main  program  and  the  ker  and  kei  functions  of  order  zero  and  one  were 
calculated  using  an  asymptotic  form  of  series.  After  the  insertion 
of  the  "BSL"  subroutine,  the  program  computed  the  Bessel  functions 
correctly . 

After  the  necessary  Bessel  functions  were  computed,  a  4  X  5 
matrix  was  set  up  involving  the  four  constants  A,  B,  C,  and  D.  Sub¬ 
routine  "SOLVE"  was  called  upon  to  solve  the  matrix.  The  last  part  of 
the  program  is  the  calculation  of  moment-curvature  relationship  from 
Eqn .  1.3-11. 

The  program  was  initially  written  in  single  precision  (i.e.  using 
eight  significant  figures)  however,  it  was  soon  discovered  that  due  to 
computer  "round-off"  misleading  results  were  obtained.  The  program  was 
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then  changed  to  double  precision  (i„e.  using  16  significant  figures) 
resulting  in  much  improved  results. 

To  facilitate  in  reading  the  computer  program,  some  of  the  follow- 


ing  symbols 

were  used: 

GAM 

- 

y,  Euler's  constant 

YM 

- 

E,  Young's  modulus 

FMU 

- 

y,  Poisson's  ratio 

TO 

- 

tQ,  amount  of  taper 

FC 

- 

c,  edge  thickness 

FK1 

- 

k,  dimensionless  parameter 

FB 

- 

b,  one-half  the  plate  width 

BETA 

6  =  c 

2t0 

FLA2 , 

FLA4, 

etc . 

- 

A2>  A*4?  etc. 

EP 

- 

e 

AL 

- 

a 

ETA 

- 

n 

BROO 

- 

ber  r\ 

BRIO 

- 

ber 4  n 

FKI11 

- 

kei£  p 

EBI02 

- 

bei2  e 

AKR01 

- 

ke^  a 

X(l), 

X(2), 

X(3) ,  X(4) 

- 

A,  B,  C,  D  respectively 

FI, 

F26 

- 

various  parts  of  Eqn.  1.3-11 

FM 

- 

moment 

R1 

— 

1/R,  curvature. 

The  input  data  consisted  of  YM,  FMU,  TO,  FC,  FK1,  and  FB  and  the 


output  was  R1  and  FM. 


. 
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C  PROGRAM  PLATE 

COMPUTER  PROGRAM  FOR  THE  EVALUATION  OF  THEORETICAL  MOMENT- 
CURVATURE  RELATIONSHIPS  FOR  PLATES  WITH  CAMBERED  B  I  TRAPE ZO I  DAL 
CROSS  SECTIONS 

COMMON  A,X,FKR00,FKI00,FKR01,FKI01,ETA,N,NI 

DOUBLE  PRECISION  A ( 4 , 5 ) , X ( 4 ) , YM, FMU* TO, FC , FK 1 , F B , P I , GAM , F J , FL A2 , 
IFLA4,FLA8,FLA,  BB, R, BETA, EP, AL, ETA, T, SUM , F I , SUMI , SUM2 , B , C , P S I , FK , D , 
IPS  I  I  ,PSI2, AI,EP2,EP3,EP4, AL2, AL 3 , AL4 , FMU2 , T03 , FM , F A , F I ,F2 ,F3,F4,F5 
l,F6,F7,F8,F9,F10,FII,FI2,F13,F14,FI5,F16,FI7,F18,FI9,F20,  F2I,F22  , 
IF23, F24,F25,F26, A23, B23,A24, B24, A25 , B25 , A26 , B26 , X 12 , X22 , X32 , X42 , 
l  BI02,A8I02,EBI02,FKI02,EKI02, AKI02, X 13 , X24, X 14 , X23 , R I , DL 
DOUBLE  PRECISION 

I  BROO,  BRIO,  BR20 ,  BROl,  BRII,  BR02,  BIOO,  BILO,  BI20,  BIOl,  Bill, 
1EBR00,EBR10,EBR20,EBR01,EBR11,EBR02, EBI 00, EB I 10 ,EBI 20, EBIOl , EBI 11 , 
IABR00,ABR10, ABR20, ABR01 , ABRil, ABR02* ABIOO, ABI 10,ABI20, ABI01 ,ABI 11 , 
1FKROO,FKR10,FKR20,FKROI,FKR11,FKR02,FKIOO,FKI 10 ,FKI 20,FKI 01 ,FKI 11 , 
lEKRO0,6KR10,EKR20,EKR01,EKRll,EKR02,EKI00,EKIi0,EKI20,EKIOl,EKIU , 
1AKR00, AKR10, AKR20, AKR01, AKR11, AKR02, AKIOO, AKI10,AKI20,  AKI01 ,AKI  11 
11  FORMAT ( 1HJ, 1F13.4»1F13.1) 

27  FORMAT (6X»6D10.3) 

C=DSQRT ( 2D0 ) 

P  1  =  3 .141 5926535 8979D0 
GAM=. 5772 1566490 l 53D0 
DC  28  J J=1 , 25 
WR  I T  E  (  6 , 29  ) 

29  FORMAT* 1HK,5X,60H  YM  MU  TO  C  K 

1  8  ) 

READ(5*27)YM*FMU,T0»FC»FK1*FB 
WRITE<6,27)YM,FMU,T0,FC,FK1,FB 
BETA=FC/(2D0*T0) 

BB=FB**2 
FMU2=FMU**2 
WR ITE ( 6 , 26 ) 

26  FORMAT < 1HL , 26H  CURVATURE  MOMENT) 

DO  28  J=l, 100,1 
F  J  =  J 

R 1=  .00 1D0*F J 
R= 1D0/R 1 

FLA4=3D0*FB»»4* ( 1D0-FMU2) / ( { T0*R )**2 ) 

FLA2=BB*DSQRT ( 3D0* ( 1D0-FMU2) )/(R*T0) 

FLA=SQRT ( FLA2 ) 

FL A8  =  FL  A4**2 
EP=2D0*FLA*DSQRT( BETA) 

AL=2D0*FLA*DS0RT( 1D0+BETA) 

ET  A=EP 
L=  1 

C  BER  ETA 

23  T=-( .2500) **2*ETA»*4/4D0 
SUM=  1D0+T 
F  I  =  IDO 
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( SA J3 ) If. 2=A J3 
5  «  *  -A  A  J  1  =  8  A  J  3 

(  a  t  j  a )  t  ><  j  a  0  *  a  j  3  *00  s = q  3 

(  ■  T  '3  +  001  )  I  >i  v, :  a  0  A  J  3  *  0  3  S  =  )  A 

q3=AT3 
1  =J 

AT3  S3H  3 

:1A\A*«.  T3*S«  »  (  )0<A.  )  -- T  t  S 

T+001 =MJ2 

001  =  1  i 


noon 
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1  FI  =  F I  +  1D0 

T=-T*.25DO**2*ETA#*4/ ( ( 2D0*F I  * ( 2D0*F  I- IDO)  )**2) 
SUM=SUM+T 

IF(DABS(T) .GE.5C-16)  GO  TO  1 
BROO=SUM 
C  BE  I  ETA 

T=„25D0*ETA**2 
SUM=T 
F 1  =  1  DO 
3  F I=F I+iDO 

T=-T*(.25DO*ETA**2)»*2/( ( < 2D0*F I - IDO ) * ( 2D0*F 1-2 DO ) )**2) 
SUM=SUM+T 

IF(DABS(T) .GE.5C-16)  GO  TO  3 
B 1 00  =  SUM 

C  BER  AND  BE  I  ETA  FOR  NU=1 

SUMl=CCOS( .75D0*PI ) 

SUM2=DS I N ( • 75DO*P  I  ) 

B= 1 • 25D0 
FI= IDO 

T=.25DO*ETA**2/ ( F I * <  F I ♦ IDO ) ) 

SUM1=SUM1+T*DC0S( B*PI) 

SUM2=SUM2+T*DSINCB*PI  ) 

5  FI=FI+1D0 
B=B+  •  5D0 

T=T*.25D0*ETA**2/(FI*(FI+1D0)  ) 

SUM1=SUM1+T*DC0S< B*PI  ) 

SUM2=SUM2+T*DSIN( B*P  I  ) 

IF( DABStT ) .GE.5C-16  )  GO  TO  5 
BR01=.5D0*ETA*SUM1 
BI01=.5D0*ETA*SUM2 
BR10  =  (BR01+BI0i  )/C 
BI10=(BI01— BROl )/C 
BR20=-BI00— BRIO/ ETA 
BI20=BR00-BI 10/ETA 
BRll=-<  BROO  +  BIOO)/C-BR01/ETA 
BI il=-( BIOO-BROO)/C-BIO 1/ETA 
BR02=-C  * ( BRO 1— B I01)/ETA-BR00 
BI02=-C*(BI01+BR01  )/ETA-BIOO 
C  KER  ETA 

I F ( ET  A. LE • 10D0 )  GO  TO  31 
CALL  BSL 
GO  TO  30 
31  1=1 

DL=DLOG( .5D0*ETA) 

SUM=— GAM 

T=— { .25D0*ETA**2 )**2/4D0 
1=1  +  1 
F  l  =  I 

PS  I =-G AM 
N=2* 1-2 
DO  14  K=1 »  N 
FK  =  K 
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oni  + 1  i-n  i 

)  #  i  *  .  <;  )  )  \  a  »  T  :i  *  *  •'  '  I  ~  =  1 

T  +  .U2=VJ2 

i  t  c :  (di-oe.  i  i.  n  )  a.HAo )  hi 

HU2  =  )Oflfl 
A  T3  138  3 

M  *  A  T  3  *  00  2  S  •  =  T 
T=MU2 
031=13 
031+13=13  l 

»  (  I  - 1  .  » OCi  Si  )  •  (  0!  i  -  1  '■*  0  )  )  )  \  X  *  a  {  , »  #  .  (  .  :  > ,  ;j  e v>  .  )  *  I  -  =  I 

7  +  mu2 -  mj2 

T  M  (dl-jc.  0.  (T  )Z ...  MI 

MU2=00IH 

i  * tw  933  Aid  I3H  31/ A  9  3H  3 

(  I  Si  5  r\  .  )  2  ;;  j  —  ,  v‘  J 2 
(  i  9  »  OCI  V  .  )\AIZ  i=SMlJ2 

oaes.i=a 

oa  i  =  n 

(  (  i Ci  1  +13  )  «  I  •  )  Y*  *  *  A  T  3  #  0  i  .=T 
ns*e  1  2  Q 0 3  #  T  ♦  I M U 2  =  1  M  32 
(  1  9  a  1  Ml  J#T  +  SMU2  =  S,MU3 

ooi  +  i  i  =  i3  e 
002. +8=8 

{  (  M  M  )  a  I  |  )  \  s  «  a  a  T  *  0  v  1 '  •  S  •  *  T  =  T 
(I  r*H  )c  1)33*  1  +  I  M2=  1MU2 
(  I  q  :  )  I  2  0*  [  +  S  12  =  Y  U  ; 

,  I  r  {  !  -  :c.  It  .  (  T  )  2 HA  ;  )  3  I 

i  2  »  A  T  1*0  G  c  .  =  i  0  ’  ) 

A  i  J 2*  A  T  3  « OQ  2  •  =  1  0 1  8 
:a {  ioi h  +  ios  0=  isfl 
;  \  ( i  - J  r  )  =  i  i 
a i  \ o i  >  •  -  i  -)-=  jS  •  ; 

i  uxoiia-  o a >■'  =  o s i a 

t  I  .  \  i  0'  00!  H  +  a  v  )  -  *i  I  -  8 

AT  :  \  I  0  1  -0  \  (  0  0  a-  I  H  )-  -  i  i  1 
J  ‘  -AT  .  \  (  i  0  I  0-  I  OVH)  *  J-  =  SOH  ; 

.il  i  -  A  I  \{  j(  <8  +  1  018  )O-  =  S0I  'J 

AT3  H3X  0 

K  OT  00  (0301  .  i  J  •  A  T  3  )  3  I 

J28  JJAO 
OE  OT  00 

i=I  I£ 

(  A  T  (»0  .2  .  )  03 J  i  =  JO 
MA0---mU2 
A  \  *  (  *  f  3  *  0  3  c  X  .  )  -  =  T 

1+1=1  ei 

1=13 
3Ac -= I  29 
S-I *S=H 
l*,i=*  Ai  03 
/i  =  X  1 
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14  PS I=PS I +1D0/FK 
D=T*  PS  I 
SUM=SUM+D 

T=— T  *( .25D0*ETA**2 )  **2/ ( { 2D0*F I* ( 2D0*F I-1D0 )  )**2) 

IF (DABS(D) .GE.5D-16)  GO  TO  15 
FKROO=— DL*BR00+.25D0*P  I*BIQO+SUM 
C  KEI  ETA 

1=1 

SUM=  {— GAM+  1D0)*.25D0*ETA**2 
T=-  ( .25D0*ETA**2)**3/36D0 
16  1=1+2 
PS I=-GAM 
DO  17  K= 1 f I 
FK  =  K 

17  PSI=PSI+1D0/FK 
D=T*PSI 

SUM=  SUM  +  D 
F 1=  I 

T=-T*(.25D0*ETA**2)**2/( (  ( F  I +2D0 ) * ( F I + IDO ) )**2) 

IF(DABS(D) .GE. 50-16)  GO  TO  16 
FKI00=-DL*BI00-.25D0*PI*BR00+SUM 
C  KER1  AND  KE 1 1  ETA 

1  =  1 

SUM1=DC0S( .75D0*PI )  * (-2D0  +  GAM+1D0 ) 

SUM2=DS I N { .7500* PI ) * (-2D0*GAM+1D0 ) 

T=ET  A**2/8D0 
B= 1. 25D0 
20  1=1+1 
F  I  =  I 
N=  I  - 1 
PS  1 1=-GAM 
DO  18  K= 1 1  N 
FK=K 

18  PSI1=PSI1+1D0/FK 
PS  1 2=-G AM 

DO  19  K= 1 ,  I 
FK=K 

19  PSI2=PSI2+1D0/FK 
PSI=PSI 1+PSI2 
A1=DSIN ( B*P I ) *PS I*T 
D=DCOS ( B*P I ) *PS I *T 
SUM1=SUM1+D 
SUM2=SUM2+A1 

T=T*.25D0*ETA**2/(FI*(FI+1D0) ) 

B=B+.5D0 

IF(DABS(D).GE.5C-16)  GO  TO  20 

FKR01=.25DO*ETA*SUMi+DCOS( . 75D0* P I ) /ETA-DL*BRO 1 +. 25D0*P I *B I 0 1 
FKI01=.25D0*ETA*SUM2-DSIN(.75D0*PI ) /ET A-DL*B 1 0 1-. 25D0*P I * BROi 
C  SOLVE  FOR  KER-KEI  PR  IME , KER-KE I  2PR I  ME , KER  1-KE I  1  PR  I  ME , KER2-KE I  2 

30  FKR10=(FKR01+FKI01)/C 
FKI10=(FKI01— FKR01)/C 
FKR20=-FK I OO-FKR 10/ETA 


i  -■  I  0  * 

i  .  ■  >'  •  * 

■  -  S. :  . 


I  29  »  T  =  J 
0  +  MU2  .12 

(5  (  (  ...  -  H  *  C  i ..  )  *  I  i » 0  wSi  1  )  \  s  *  *  (  5  3  A  I  *  ■  ’ G  5  .  )  *  1  -  T 

cl  T  )0  {3J-32.J0.  (0)2  A  1)  I 
-*  2-1  JO  '  IdS.+O  .  ■  .  •  J  J  -  =  0  ■’  X  • 

A  T3  13X  3 

1  =  1 

■  *«AT  i»  03  ..a  (0:1  ♦  '*1 A  <)  ~  )  =  ^  U  2 

0(]  )t  \i  **  (  5*»AT  «  00  2  ‘  .  )  -  =  T 

5+1=1  31 
MAO-= 1 29 

i  *  i  =x  vi  oa 

X  =  X9 

X9\00  (+129=129  VI 
I  29  «  T  =  ! 

C1  +  MU2=MU2 
1*13 

(  *  *  {  {  I  +  1  .  )  *  (0  >  +  I  i  )  >  )  \  f .  *  *  (  S  a  a  A  ■;  *  0  ies  .  )*  .'-  =  *! 

31  T  )0  (  31  -  Jc  .  .(D2  AQ  )  -i  I 

mU,  +  •  I  <  »  .  .  Ida  J  I--  AO  I  1 

AT  3  I  I  3  >1  QUA  j  83X  0 

i  =  I 

(0-, .[  *  a: )V- )  » (  1  9  *002N  .  )  o  j=  l  MJ2 
{  0  1  +  1 5 -  )  *  (  l  9  *  Ik'V  .1/12  J 2 

G'AS*»A  |  j-T 
0025  . 1 =a 

1+1=1  05 
1  =  19 
1-1=/ 

• AO-=i 1 29 

/ , ; =  x  8i  on 

X  =  X9 

X9\03J +1 129=1 129  dj 
AQ--  S 1 29 
I  » 1  =  X  PI  00 
X  =  X9 

X9\0 01+51 29=5129  Pi 
5  I  29+ I  129  =  129 
T  » 1 1  9 *  (  19+  )  \A\  2  0  =  1  A 

T*  1  29* (  ! 9*8 ) 23 JO  =  0 

a+ lMU2=i*U2 

I A+SMU2=SMU2 
(  (  .  ■ ;  1  +  I  9  )  *  ]  9  )  \  5  ’  »  a  T  ) « i  :  2  s .  *  T  =  1 

002. +8=9 

05  »?  On  ( 3l-n  . JO. (0 )28A j  )  91 

■i .  ♦  10  h >  i  -■  Ki-  AT  \  (  1  9  » 0 G c  V  .  )  2 U 00  +  I  -  U  <  *  A  T  J  *  )0 <  5  .  =  f  •  X9 
S  *  —  1 0  1  6  *  .1  • !  -  A  1  .:  \  (  i  '■.*(■:!  ;.M  .  >1-^1  -  »2*AT3«  ••101 

,  :  ]  a-i  ,  :  .  i  I  *0H  3 V J  3 2  0 

j\ (  !  0  1  X  9+1 08X9 ) =0] V*  1  Ot 
0\ ( 1 09X9-1  )I  ;9  1=011X9 
A  r  ■  V  i  9X  i  -00  I  X  9  -  =  )5>|X9 
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25 


FKI20=FKR00-FKI 10/ETA 

FKR1 l=-(FKR00+FKI00 )/C-FKR01/ETA 

FK  1 1 1=— ( FK  I 00— FKR00 )/C-FKI01/ETA 

FKR02=-C*{FKR01-FKI01)/ETA-FKR00 

FK  I  02=-C*  (FKIOl  +  FKROD/ETA-FKIOO 

I F ( L  .EQ • 2 ) GO  TO  25 

EBR00=8R00 

EBR10XBR10 

EBR20=BR20 

EBRO  1=BR0 1 

EBR11=BR11 

EBR02-BR02 

EB 1 00=B 1 00 

EB 1 10  =  B  I  10 

EB  1 20=B I  20 

EB 1 0 1  =  B 101 

EB  1 1 1=B  I  1 1 

EB  I02-B 102 

EKROO=FKROO 

EKR 1 0=FKR 10 

EKR20=FKR20 

EKRO  1=FKR0 1 

EKR 1 1=FKR 1 1 

EKR02-FKR02 

EK I 00=FK 100 

EKI10=FKI10 

EK I 20-FK I 20 

EKI01=FKI01 

EKI11=FKI11 

EK  I02=FK 102 

ETA=AL 

L=2 

GO  TO  23 
ABROO=BROO 
ABR 1 0=BR 10 
ABR20=BR20 
ABRO 1=BR0 1 
ABR 1 1=BR 1 1 
ABR02=BR02 
AB I 00=B 100 
ABI 10=BI 10 
AB 1 20=8  1 20 
AB I01=B 101 
ABI il=BIll 
AB 1 02=B 1 02 
AKROO=FKROO 
AKR 1 0=FKR  10 
AKR20=FKR20 
AKRO 1=FKR0 1 
AKR 1 1=FKR  1 1 
AKR02=FKR02 
AK I 00=FK 100 


I  3  \  .  I  I  X3-  )'■  1  =  .;  i  I  ■:  ■ 

!  \ .  )  -  \ (  oi  1  y  ■  • ; -  ■  i 

r  d \  i -  (  I-  f  ■ ;  -  ■ 1  .  i 

■  -AT3VU  I  -  :  }#=U*D-*S 

I  -  I  \  (  1  -  H  ,'•)••  - 

es  UT  OOCS.tH.  J)3I 
00X0=00 

oi xa=o i aoa 

0$*8=OSXH3 
i 0M8  =  i 0X03 

i  i«a=i  ifiir) 
so.<  a=soxa3 
ooi a =oo i  h 

01 1  a =01 1  a  3 
os ia=osi a a 
i  o  i  a  =  i  o  p  3 
ii 18=1 i IH3 
S0I8=£01 8 
00XX3=Q0XX3 
0iXX3  =  :jiXX3 
osax  i=oxx  j  1 

I  0MX3=1  OHXfi 
i  ix*3=i ixxa 
S0MX3=S0XX3 
001  X  3=00  I  *  3 
01  I X3=0i I x  « 
osix3=osi xa 
i  o  i  x  a  =  i  u  l  x 
I  1  1X3  =  1  1  I  * 3 
S0IX3=S')I  X3 
J  A  =  A  f  3 
$  =  J 
£S  JT  00 

ooxa=ooxaA  es 
oixa=oixaA 
OSX0=OSX0A 
ioaa=i omha 

I  1  '•  H=i  1X4A 
S0X3=S0flHA 

ooi a=ooiflA 
oi  i a = o 1 1 ha 
OS  I H  =  OS I  HA 

ioia=i uiha 

1 i I H  =  i  i  I  HA 
SOI «' i  —  S  )I8A 
00XX3=00HXA 

oixxh=oihxa 

0\xx3=0<P  xa 
1  0>  X1--I  0«  XA 
1  1  fl‘  31  1 HXA 
SOX  •  3  =  S 0 M X A 
00  I  X 3  =  0  01  x/ 
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AKI 10=FKI  10 

AK 1 20=FK 1 20 

AKI0i=FKI01 

AKI11=FKI11 

AKI02=FKI02 

EP2=EP**2 

EP3=EP**3 

EP4=EP**4 

AL2= AL**2 

AL3=AL**3 

AL4=  AL  *  *4 

At  1,  1)=-EP3*EBR00-24D0*EP*EBI00-48D0*EBR10+8D0*EP2*EBI 10 
A ( 1,2)=-EP3*EBI00+24D0*EP*EBR00-48D0*EBI 10-8D0*EP2*EBR10 
At  1*  3)=-EP3*EKROO-24DO*EP*EKIOO-48DO*EKR10+8DO*EP2*EKI 10 
At 1,4)=-EP3*EKIOO+24DO*EP*EKROO-48DO*EKI 10-8D0*EP2*EKR10 
At  1 ,5)=0D0 

At  2, 1)=-EP2*EBI10+4D0*EP*EB100+8D0*EBR10 
A (2,2)=  EP2*EBR10-4D0*EP*EBR00  +  8D0*EBI  10 
A(2,3)=-EP2*EKI 10+4D0*EP*EK I 00+8D0*EKR 1 0 
At  2,4)=  EP2*EKR10-4D0*EP*EKR00+8D0*EKI10 
A(2,5)=FMU*BB*EP3*BETA/(R*FLA2) 

At  3, 1)=-2D0*ABR10-AL*ABI00 
A(3,2)=-2D0*ABI 10+AL*ABRO0 
A(3,3)=-2D0*AKR10-AL*AKI00 
A(3,4)=-2D0*AKI 10+AL*AKR00 
A ( 3 , 5 )  =-T  0  * (  1D0+FK1) *AL3/( 2D0*FLA2) 

A ( 4 , i)=AL2*ABR00-2D0*AL*ABI 10-EP 2*EBR 00+2D0*EP* EB I  10 
A(4,2)=AL2*ABI00+2D0*AL*ABR10-EP2*EBI00-2D0*EP*EBR10 
A(4,3)=AL2*AKR00-200*AL*AKI 10-EP2*EKROO+2DO*EP*EKI10 
A(4»4)=AL2»AKIOO+2DO*AL*AKRiO-EP2*EKlOO-2DO*EP*EKR10 
A(4, 5)=— 6D0*FMU*BB/ ( 4D0*R *FL A4 ) * ( AL4-EP4 ) 

N=4 

CALL  SOLVE 
T03=T0**3 

FM=YM*FB*T03/ ( ( 1D0-FMU2 ) *768D0*FL A8*R  ) * ( EP**8-AL**8) 
FA=YP*FMU*T03/(  t 1 2D0 *FB*FL A4 ) * ( 1 D0-FMU2 ) ) 

F 1 = ( EP2  *EBR00— 2  DO*EP*EB 1 10— AL2*ABR00+2D0*AL*ABI 10)*X( 1 ) 
F2=(EP2*£BI00+2D0*EP*EBR10-AL2*AB I 00— 2D0* AL* ABR 10 ) *X ( 2 ) 
F3= ( EP2»EKR00-2D0*EP*EKI10-AL2»AKR0Q+2D0*AL*AKl 10) *X( 3) 
F4=( EP2*EKI00+2D0*EP*EKR10-AL2*AKI00-2D0*AL*AKR  10 )*X(4) 
F5=-(EP3*EBR10-AL3*ABR10 )*X( 1 )+3D0*Fi 
F6=-(EP3*EBI 10-AL3*ABI10)*X( 2)+3D0*F2 
F7=-(EP3*EKR10-AL3*AKR10)*X( 3)+3D0*F3 
F8=- (EP3*EKI 10-AL3*AKI 10) *X( 4)+3D0*F4 
FM=FF+FA*3D0/2D0*(F1+F2+F3+F4+F5+F6+F7+F8) 

F  9=(EP4*EBR20-AL4*ABR20)*X( 1)+4D0*F5 
F10= (EP4*EBI20-AL4*ABI20) *X( 2 ) +4D0*F  6 
FI  1= (EP4*EKR20-AL4»AKR20) *X( 3)+4D0*F7 
F12= (EP4*EKI20-AL4*AKI20)*X(4)+4D0*F8 
FM=FM+FA*.5D0*( F9+F10+F11+F12) 

FM=FN+9D0*YM*FMU2*FB*«5*T0/( 2D0*R**3*FL A2**6 ) *  t EP4-AL4) 
F13=F1 
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S 0 l X3 - S 0  I  X  A 
S*#93=S93 
t  *  * ( !  i~  t  S  : 
A  »  *  M  J  =  A  M  -I 
S*  *  JA  =S  IA 

J  >  *  JA  r.  f  JA 
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1  •••'.*  j /■  -01  ft  ./  *(Jl\-=(  i  «  v  )  A 

00.  «  !  a ♦ o i  1  I  A »  -=(S,£  )A 

01  XA»  JA-  J  1A  »  )S-*(f  .  e  ) 

0  "  (  )■  *  .  :>  f  O' I  I  »0.’S-=  (  .  t  •  ) 
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»  I  *<  .  *  ■  »S93-(  4  i  HA*  JA*0GS-009XA*S  I  *  (I  *A) 
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F14=F2 
F15=F3 
F 16=F4 

FM=FM-6D0*YN*FMU*FB**3*T0/( R«*2*FLA8 )*(F13+F14+F15+F16) 

X 1 2=X { 1 )**2 
X22=X ( 2 ) **2 

F17= (X12  +  X2  2  )*( EP*( EBRO  1* FBI  1 1~E BR 1 1 *EB I Oi )-AL* ( ABRO 1 » AB I  1 1-ABR 1 1  * 
1ABI01) ) 

F18= (X12-X22 } *< EP2*( 2D0*E BRO 1*EB 1 0 1-EBROO*EB I 02-EBR02*E B I  00 )- AL2* { 
12D0*ABROi*ABI01-ABR00*ABI02-ABR02*AB 100) ) /200 
X32=X ( 3 ) **2 
X42-X ( 4  )  **2 

F19=(X32+X42) *( EP* < EKRO 1*EK I  1 1-EKR11*EKI 01 )-AL * ( AKRO 1  * AK 1 1 1-AKR 1 1 « 
1AKI01) ) 

F2O={X32-X42)*(EP2*(2D0*EKR01*EKIOl-£KROO*£KlO2-EKR02*EKI00)-AL2* ( 
12DO*AKROl*AKI01-AKR00*AKI02-AKRO2*AK 100  )  )  /2D0 
F21=-X( 1 )*X (2  )*(EP2*(EBR01**2-EBR00*EBR02-EBI01**2  +  EBI00*EBI02 )-AL 
1 2* ( ABRO 1 **2~ ABROO«  A8R02-AB I  01**2 +A BIO 0*ABI 02) ) 
F22=-X(3)*X(4)*(EP2*{EKROl**2-EKR0O*EKR02-EKIOl**2+EKlOO*EKIO2)-AL 
12*(AKR01**2-AKR00*AKR02-AKI01**2+AKI00*AKI02) ) 
A23=EP2*(2DO*EBROl*EKIOl“EBR0O*EKIO2-EBRO2*EKI00+2DO*EBI01*EKROl-E 
1BIOO*EKR02-EBI02*EKROO) 

B23=AL2*(2D0*ABR0i*AKI01-ABR00*AKI02-ABR02*AKI00+2D0*ABI01*AKR01-A 

1BI00*AKR02-ABI02*AKR00) 

X 13=X ( 1 )*X(3) 

X24=  X  (  2  ) *X ( 4  ) 

F23=(X13-X24 )/2C0*( A23-B23) 

A24=-EP*(EBRll*EKI0l-EBR01*EKIll-EBI 1 1*EKR0H-EB  I01*EKR  1 1 ) 

B24=-AL*<  ABRli*AKI01-ABR01*AKI ll-ABI l 1* AKR01+ AB I 0 1* AKR 1 1 ) 
F24={X13+X24)*( A24-B24) 

X 14=  X ( 1 ) *X ( 4  ) 

X2  3=  X ( 2 ) *X (  3  ) 

A25=EP2*{-2D0*EBR01*EKR0H-EBR00*EKR02+EBR02*EKRG0+2D0*EBI01*EKI01- 

1EBI00*EKI02-EBI02*EKI00) 

825=AL2*(-2CO*ABR01*AKROl+ABROO*AKRO2+A6RO2*AKR0O+2D0*ABIOl*AKIOl- 

1ABIOO*AKI02-ABI02*AKIOO) 

F25= (X14  +  X2  3 )/2C0*<  A25-B25) 

A26  =  EP*(EBR11*EKR01-EBR01*EKR1H-EBU 1 *EK 1 0 1-EB I 01*EK I  11) 
B26=AL*(ABRli*AKR01-ABR01*AKRll«-ABIll*AKI0l-ABI01*AKIll) 

F26= ( X14-X23 ) * ( A26-B26 ) 

FM=FM+YM*FB*T0/(4D0*R*FLA4 ) * ( F 17  +  F 18 +  F 19+F20+F2 1+F22+F23+F24+F25+F 
126) 

28  WR I TE( 6  « 1 1  )  R1,  FM 
STOP 
END 
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SUBROUTINE  SOLVE 

COMMON  A,X,FKR00,FKI00,FKR01,FKI01,ETA,N,NI 

OOUBLE  PRECISION  A { 4 , 5 ) , X < 4 ) , FKR GO , FK I  00 , FKROi , FK I  01 , E TA , AM AX , TS , 
IFMULTPtFNUM 
M=N+  1 
NM1=N-1 
00  100  1  =  1  ♦  N 

100  A(I,M)=-A(I,M) 

DO  101  1=1, NM1 
AM AX= A (  1,1) 

K=  I 

DO  102  J=  I  ,  N 

IF(DABS(A( J,I )  )  .LE.DABSI AMAX) )  GO  TO  102 
AM AX=A ( J , I ) 

K=J 

102  CONTINUE 

I F I AMAX.NE • 0 * 0 ) GO  TO  103 
RETURN 

103  IFIDABSI AMAX ) .GE.  ID-21 )  GO  TO  104 
WRITEI6, 105) 

105  FORMAT ( 1HJ , 14X  *  45H  MATRIX  IS  EITHER  SINGULAR  OR  ILL  CONDITIONED) 

104  DO  106  J=I,M 
TS= A (I , J ) 

A ( I  ,  J)  =  A  ( K »  J  ) 

106  A<K,J)=TS 
IP1=I+1 

DO  107  J= I P 1 ,  N 
FMULTP= A ( J , I  l/AII  ,  I  ) 

DO  107  L= I P 1 »  M 

107  A( J,L)=A( J,L )-FMULTP*A( I,L  ) 

101  CONTINUE 
X(N)=A(N,M)/A(N,N) 

DO  108  1=1, NM1 
K=N-  I 

KP 1=K+ 1 
FNUM=ODO 
DO  109  J=KP 1 ,  N 
109  FNUM=FNUM+X(J)*A(K,J) 

108  X(K)={ A(K,M)-FNUM)/A(K,K) 

RETURN 

END 


V  J  U  C  1  T  !  J  U  G 


1  ,  .  ,  r  ;i  t  i  0  I  >  i  t  U.I  >H  ♦  G  I  OV  i  t  •<  t  A  /  '  40 
t  t  f  ;  ,!  1  ,  ,  1  ,1,:.;  1  >1,00H>  *  (  A  )  v  ,  ( <  ,  t  )  '  !  *.' 

MUH9,9TjUM  U 
1  +H*M 
I-H  =  i AA 
Af 1=1  001  00 
(  t 1 ) A~=(  « I ) A  001 

1  Ml/I ,  1  =  1  101  00 

(1,1)  -  ■ A  M  A 

A  I  =  X 

H ,  I  =1  SOI  )0 

A  I  i  r  10  U  X  A  A  )  t  A  .•)..((  1  t  L  )  A  )  (i  HI 

(1,1  }  A  =  X  A  M  A 
l=X 

-tUHI  TH03  SOI 

>■  1  or  f:  u  (0.0.  90  .  X  ^  M  A  H  I 

/i«UTd« 

-01  OT  :  (  IS-Ol.jO.  (  X  A  M  A  )  £  8  A  I  )  ~U  £01 

(501 td)3TIfltf 

■  .  ...  j  .  ;  I  .  ■  I  I  r  AM  H  tXAltlMJ  >?J  i  I 

M,I=L  001  10  AO  1 

( L , I )  A  =  2  T 
(  l ,  >i )  A  =  (L  ,  I  )  A 

2 T  = ( t , »  1  A  dOl 
1+1=191 
H  ,  i  9  I  =  L  \0 1  00 
(  I  ,  I  )  A  \  (  1  ,  L  )  A  =9  T  JU  M 
M , IS  I =J  TO 1  00 

(  J  ,  I  )  A  *  9  T  JUH-  (  J  ,  U  A  =  (  J  ,  L  )  A  Hi 

3UHITH33  1  l 
( H  ,  A  )  A  \  (  M  , ,  )  A  =  (  0  )  x 
i  MU,  1  =  1  .101  00 

I  -/=* 

1 +X= 1 9M 
000= MU HO 
H,  19  >1  =  1  001  30 
(  t.  ,  x  ) ;  *  (  L  )  X  +  UO  1=  J.  9  GO  i 
( X ,  ) A\{  ‘  J :  H  - (  *i  X  )  A  )  =  (X ) X  b  01 

/ HU T  d 9 
0H3 


o  o 


SUBROUTINE  BSL 


ASSYMPT  OTIC  EXPANSION  OF  KER  KEI  KER I  KEI1 
COMMON  AtX»FKR00»FKI00fFKR01fFKI0I*ETA»N*NI 

DOUBLE  PRECISION  A ( 4 , 5 ) , X ( 4 ) , FKROO ♦ FK 1 00 , FKRO 1 , FK I  01 , E T A , S2 , SUM  1 , S 
1UM2,PI ,BETA,FC,C,S,B,T,FI 
S2=DSQRT(2D0 ) 

P 1=3 *1415926535897900 

SUM  1  =  IDO 

SUM2=000 

8ETA=ETA/S2+PI/8D0 

FC=DSQRT(PI/(2D0*ETA )) *DEXP < -ETA/S2 ) 

B=  *2500 
1=1 
JC  =  2 
JS  =  4 
JC  P= 4 
JSP=2 
JCCP=8 
JSSP=6 

T=-1D0/<8D0*ETA) 

SUM1=SUM1+T*CCQS(PI*B) 

SUM2=SUM2+T*DSIN(PI*B) 

3  1=1+1 
F I  =  I 

B=B+  *25D0 

I  F (  I  *EQ.  JC  )  GO  TO  4 
I F (  I  *EQ* JCP  )  GO  TO  22 
IF (  I  *EQ*  JCCP )  GO  TO  23 
C=DCCS ( P I *B ) 

GO  TO  5 

4  C=0D0 
JC=  JC+4 

T=T*(-(2D0*FI-1D0)**2)/(FI*8D0*ETA) 

GO  TO  18 

22  C=-100 
JCP= JCP+8 
GO  TO  5 

23  C=  IDO 
JCCP= JCCP+8 

5  T=T*(-(2D0*FI-1C0)**2)/(FI*8D0*ETA) 

SUM1 =SUM 1+T*C 

18  IF(I.EQ.JS)  GO  TO  6 
I F (  I  *EQ  .  JSP  )  GO  TO  24 
IF ( I*EQ.JSSP)  GO  TO  25 
S=DSIN( PI*B) 

GO  TO  7 

6  S=ODO 
JS= JS+4 
GO  TO  19 
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24  S= IDO 
JSP= JSP+8 
GO  TO  7 

25  S=-1D0 
JSSP=JSSP+8 

7  SUM2=SUM2+T*S 

19  IF(DABS(T) .GE.1D-08)  GO  TO  3 
C=DCCS ( BET  A  ) 

S=DS  IN { BET  A  > 

FKROO=FC*(SUMl*C— SUM2*S) 
FKI00=FC*(-SUM1*S-SUM2*C) 

1=1 
JC  =  2 
JS=4 
JCP=4 
JSP=  2 
JCCP=8 
JSSP=6 
SUM 1= IDO 
SUM2=0D0 
B=  • 25D0 

T=3D0/(8D0*ETA) 

BETA=6TA/S2+.625D0*PI 
SUMi=SUMl+T*DCOSI PI*B ) 
SUM2=SUM2+T*DSIN(PI*B) 

12  1=1+1 
F  I  =  I 

B=B+.25D0 

IF(I.EQ.JC)  GO  TO  8 
IFU.EQ.JCP)  GO  TO  26 
IF(  I.EQ.JCCP)  GO  TO  27 
C=DCCS{PI*B) 

GO  TO  9 

8  C=0D0 
JC= JC+4 

T=T*(4D0-(2D0*FI-lD0)**2)/( FI*8D0*ETA) 
GO  TO  20 

26  C=-1C0 
JCP= JCP+8 
GO  TO  9 

27  C= IDO 
JCCP= JCCP+8 

9  T=T* (4D0-{2D0*FI-1D0>**2) /( FI*8D0*ETA ) 
SUM 1  =  SUM 1  +  T  »C 

20  IF(I.EQ.JS)  GO  TO  10 

I F (  I.EQ.JSP)  GO  TO  28 
IF ( I *EQ* JSSP )  GO  TO  29 
S=DSIN(PI*B ) 

GO  TO  11 
10  S=0D0 
JS= J  S  +  4 
GO  TO  21 
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28  S=1D0 
JSP= JSP+8 
GG  TO  11 

29  S=-1C0 
JSSP= JSSP+8 

11  SUM2=SUM2+T*S 

21  IF(OABSIT).GE.1D-08)  GO  TO  12 
C=DCOS( BETA) 

S=DS  IN ( BET  A ) 

FKR01=FC*(SUM1*C-SUM2*S) 
FKI01=FC*(-SUM1*S-SUM2*C ) 
RETURN 
END 
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APPENDIX..  B 

COMPUTER  PROGRAM  FOR  EVALUATING  BENDING  MOMENT.  AND  LONGITUDINAL  CURVATURE 


FROM  STRAIN  GAUGE  READINGS 


The  program  shown  below  determined  the  curvature  using  Eqn.  2.4-2. 
The  particular  program  shown  is  for  Plate  II-C.  The  program  was 
applicable  to  all  the  plates  tested  after  the  appropriate  changes  were 
made  to  it.  The  input  data  used  was  the  double  bending  strain  e^. 

After  the  radius  of  curvature,  R,  was  calculated,  the  program  determined 
the  correction  moment,  M^ ,  using  Eqn.  3.2-1.  The  correction  moment  was 
then  subtracted  from  the  applied  moment,  yielding  the  correct  value  for 
the  moment.  Output  consists  of  e  and  the  corresponding  support 
distance,  applied  load,  arc  length  from  mid-line  to  measuring  station 
under  consideration,  correction  moment,  curvature  and  moment.  For  ease 
in  reading  the  program  the  following  symbols  were  used: 


w 

-  weight  per  unit  length  of  plate 

SUPD  or  A 

-  one-half  the  distance  between  the  end  supports 
of  the  plate 

THNS 

-  distance  between  strain  gauges  at  a  particular 
measuring  station  (i.e.  thickness  of  plate  at 
centerline) 

S 

-  distance  from  mid-line  to  measuring  station 
measured  along  the  plate 

FL 

-  applied  load  to  pulley  rims 

E  (L)  or  ET  -  EB  -  e„ 

D 


FM 

-  moment 

Y 

-  P,  horizontal  distance  from  mid-line  to  measuring 
station 

FM1  or  Ml 

-  M^ ,  correction  moment 

Rl 

-  1/R,  curvature 

The  symbols  not  defined  above  are  self-explanatory  in  the  program. 
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COMPUTER  PRGGRAM  FOR  THE  EVALUT I  ON  OF  EXPERIMENTAL  RESULTS  FOR 
PLATE  II-C 


W=- 1333 

SUPD=16.75 

THNS=.221 

1  FORMAT! IX, 1415) 

2  FORMAT! 1X,7F10.3) 

3  FORMAT! 1H J , 4F 10 . 2 , 2 F 10 . 4, IF  10.  1) 

DIMENSION  J ( 340 ) , E ( 340 ) 

DO  6  1=1,5 
RE AD  !  5 , 2 )  S 
WR I TE ( 6  »  5 ) 

5  FORMAT! 1HK,70H  ET-EB  SUP  DIST  LOAD  ARC  Ml 

1  l/R  MOMENT  ) 

RE AD ( 5 , 1  )  ( J ! K ) ,  K  =  20,340,20) 

DO  6  L=20, 340,20 
FL  =  L 

E  (  L  )  =  J  !  L  ) 

FM=7.625*FL 
R=THNS/ ( E ( L ) *l.E-06) 

TH-S/ ABS ( R ) 

Y=ABS!R)*SIN(TH) 

A=ABS(R)*SIN(SUPD/ABS!R)  ) 

FMl=-W*SUPD*l A-Y)-W*ABS(R )*! SQRT { R*R- A* A ) -SORT ( R*R-Y*Y) +  Y*{ SUPD/AB 
1S(R)-TH) ) 

FM=FM+FM 1 
Rl=l./R 

6  WR I TE ! 6 , 3 )  E(L),A,FL,S,FM1,R1,FM 
STOP 

END 
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